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MINIMAL SURFACES IN THE 3-SPHERE BY STACKING CLIFFORD TORI 


DAVID WIYGUL 


Abstract. Extending work of Kapouleas and Yang, we construct by gluing methods sequences of 
closed minimal surfaces embedded in the round 3-sphere. Each surface resembles multiple tori, all 
coaxial with a nearby Clifford torus, that have been connected by many small catenoidal tunnels. 
The sequences can be indexed by the number N of tori joined and by a rectangular lattice on 
the corresponding Clifford torus, aligned with its axis circles. The terms of a given sequence 
themselves are then indexed by the sublattices of this lattice, which determine the locations of 
the catenoidal tunnels. The elements of 0(4) preserving a surface in the sequence constitute a 
discrete group preserving the corresponding sublattice, and this symmetry group contains isometries 
exchanging the two sides of the associated Clifford torus precisely when the lattice is square. A 
given sequence converges to the Clifford torus counted with multiplicity N, while an appropriately 
blown-up sequence of tunnels at a fixed location converges to a standard Euclidean catenoid. The 
finite number of tori and rational ratio of lattice spacings may be prescribed at will, but the number 
of tunnels, or equivalently the surface’s genus, must be taken large in terms of this data in order 
to complete the construction. 


1. Introduction 

In Kapouleas and Yang constructed a sequence of embedded minimal surfaces in the round 3- 
sphere converging to the Clifford torus T counted with multiplicity 2; each surface approximates 
two parallel tori coaxial with a nearby Clifford torus (not included in the surface) in between them 
and joined by many small catenoidal tunnels centered on a square lattice on the central Clifford 
torus. Accordingly they called their surfaces doublings of the Clifford torus. Kapouleas announced 
these in as the first examples of a general class of gluing constructions to double given minimal 
surfaces, subsequently discussed further in [^, and he has recently produced embedded doublings 
of the equatorial 2-sphere in 

The surfaces of Kapouleas and Yang are highly symmetric, admitting many horizontal symme¬ 
tries, which preserve as sets the two sides of the doubled Clifford torus and permute the lattice sites, 
as well as vertical symmetries, each of which exchanges the two sides of the doubled torus but fixes 
as a set a catenoidal tunnel. In fact all these symmetries are enforced throughout the construction 
and exploited to simplify its execution. The present article carries out less vertically symmetric 
doublings of the torus, with the symmetry broken in two ways. First, we allow the catenoidal 
tunnels to be arranged on rectangular rather than strictly square lattices. Any isometry of 
exchanging the two sides of the doubled Clifford torus will fail to preserve such a lattice, unless it is 
square. Second, we interpret doubling in a generalized sense, realizing also triplings, quadruplings, 
and in fact embedded minimal surfaces resembling an arbitrary finite number of tori coaxial with a 
nearby Clifford torus and connected to one another by many small catenoidal tunnels. Whenever 
at least three tori are incorporated, even if these tunnels are centered on square lattices, there will 
be many tunnels admitting no symmetries that preserve each as a set but exchange its ends, nor 
will the symmetry group act transitively on the collection of tori. 

These new constructions add to the list of known closed minimal embeddings in S^, so far 
comprising those found in d, 0, |, d, and [^. (Actually in this last reference Marques 
and Neves apply min-max methods to prove the existence of infinitely many minimally embedded 
hypersurfaces in every closed manifold having positive Ricci curvature and dimension between 3 

1 


and 7.) A variety of further examples, to be obtained by sweepouts, have been described by Pitts 
and Rubinstein [^, some of which appear consistent with surfaces obtained elsewhere (among the 
above references) by alternative techniques (but none look like the surfaces we introduce here). The 
survey article by Brendle contains an outline of several of the constructions just mentioned. 

The constructions at hand should be of interest not only as providing new examples of minimal 
surfaces in but also as a basis for further doublings with asymmetric sides. While this work con¬ 
fronts configurations with less vertical symmetry than is present in [^, Kapouleas’ sphere doubling 
addresses a situation with less horizontal symmetry. We anticipate applications demanding a 
synthesis of the two refinements. A program toward doubling constructions of increasing generality, 
including potential applications, is described in [^. 

The present work naturally emulates, with a few departures, the approach of and draws 
extensively from the general gluing technology developed by Kapouleas, much of which can be found 
summarized in and was itself inspired by techniques from 12 . Although the current article can 
be read without reference to or any other gluing constructions, for the rest of this introduction 
we will make use, without detailed explanation, of terminology standardized by Kapouleas, so that 
the reader already acquainted with it may easily appreciate the principal differences between this 
construction and [^. 

We now outline our procedure in very rough terms. As basic data we take positive integers 
N > 2, k < i, and m, where N specifies the number of copies of T to be made, while k, i, and m 
determine — 1 rectangular lattices, with ratio kji of fundamental edges, on which are centered 
in total kim?{N — 1) catenoidal tunnels. Corresponding to a choice of such data, an initial surface 
is built as follows. We start with N tori coaxial with T that have small constant mean curvatnre 
and so lie near T. From each of the two outermost tori we excise kim? discs centered on certain 
rectangnlar lattices, and from each of the remaining intermediate tori we excise 2k£m? discs. Using 
local coordinates for adapted to the tori, we take standard catenoids, truncated at a pair of 
large, opposite circles of revolution, shrink them, and glue them to the tori, along boundary circles 
of the components so as to bridge adjacent tori, bending the ends as necessary so that the resulting 
connected surface is smooth. 

For m large in terms of the other data, the tori approach T and the catenoidal tunnels, after 
rescaling, converge to a standard catenoid in Euclidean space. Thus for large m each such initial 
surface is, in a certain sense, approximately minimal, and the construction is completed by per¬ 
turbing the surface to exact minimality. Two mechanisms of perturbation are applied in tandem. 
One sort of perturbation is realized by considering graphs of small functions over the initial surface. 
To select the right function is then to solve the quasilinear elliptic differential eqnation prescribing 
zero mean curvature for the corresponding graph. The equation can be stndied by comparing the 
linearization of the operator governing the mean curvature of graphs to certain large-m limit oper¬ 
ators on the limit catenoids and limit torus. In the simplest scenario one conld solve the linearized 
operator on the toral and catenoidal components separately, combine these solutions through an 
iterative procedure, and finally invoke an inverse function theorem to solve the original nonlinear 
equation. However, the presence of nontrivial kernel to the limit operators gives rise to approximate 
kernel that thwarts the approach just described. 

The space of admissible perturbing functions is constrained to respect the symmetries enjoyed 
by the initial surface, and so their imposition has the effect of reducing the dimension of the 
approximate kernel. Each torus turns out to carry one-dimensional approximate kernel of its own, 
but in the two tori can be exchanged by reflections through certain great circles, and so together 
the tori contribute just one dimension to the approximate kernel in versus N dimensions more 
generally. Furthermore, in these reflections through circles render trivial the approximate kernel 
on the catenoidal tunnels. Following the approach of in the absence of these symmetries, each 
tnnnel wonld feature one-dimensional approximate kernel, but we bypass this kernel altogether by 
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altering, as compared to [^, the initial data at the tunnel’s waist for the rotationally invariant 
mode of the solution. 

To overcome the obstruction posed by the approximate kernel, introduces substitute kernel, 
spanned by a single function supported on the tori away from the circles where they attach to the 
tunnels. By adding multiples of this function to the source term of the linearized equation, the 
so modified source can be made orthogonal to the approximate kernel, enabling the success of the 
above scheme, but at the cost of solving the original equation only modulo substitute kernel. For 
the same purpose the current construction introduces iV-dimensional substitute kernel, spanned by 
functions each of which is supported on a single torus away from the tunnels. (Actually, in this 
construction we never explicitly identify the approximate kernel, nor do we invoke the h metric 
employed in for its analysis, but our application of substitute kernel is morally identical.) 

A further difficulty concerns the vast disparity in scale between the waist radii of the catenoidal 
tunnels on the one hand and the much greater spacing between the tori on the other. The initial 
surface’s second fundamental form grows toward the waists of the catenoids from a small value 
bounded uniformly in m on the tori to a value diverging with m, and the embeddedness of graphical 
perturbations is most precarious near the waists. For these reasons as well as to ensure convergence 
of the iteratively defined global solution, it is necessary to arrange for solutions on the tori to decay 
toward the catenoidal waists. 

All of the catenoids attaching to each of the two outermost tori—the only type of torus ap¬ 
pearing in —are equivalent modulo the symmetries, and adjustment of the source term by the 
substitute kernel suffices to achieve such decay on these catenoids. (Again, our actual approach 
deviates somewhat from this description, applicable to [^, but just superficially.) However, each 
of the intermediate tori, — 2 in number, attaches to catenoids of precisely two classes, and so 
the appropriate decay of solutions requires the introduction of another N — 2 functions, linear 
combinations of which are added to the source term to arrange decay, a device originating in 
but unneeded in [^. In total we arrive at a {2N — 2)-dimensional extended substitute kernel, the 
sum of the substitute kernel and the span of these functions, modulo which subspace we can, for 
large m, invert the linearized operator. 

Thus an infinite-dimensional problem is reduced to a finite-dimensional one. The resolution of 
this latter problem requires the second type of perturbation and is best understood in terms of a 
correspondence, which Kapouleas calls the geometric principle, between the initial geometry and 
the analytic obstructions that the extended substitute kernel represents. Coarsely put, elements of 
the extended substitute kernel can be generated, as components of the initial surface’s mean cur¬ 
vature, by certain motions of its building blocks—here catenoids and tori—relative to one another. 
So motivated, the other type of perturbation is realized by incorporating parameters, one for each 
dimension of extended substitute kernel, into the definition of the initial surface whose variation 
repositions the component tori and catenoids. Thus for each choice of k, i, m, and N, we define 
not one initial surface but a {2N — 2)-parameter family of them. 

More specifically, two parameters may be associated with each of the N — 1 classes of catenoids, 
one for each pair of adjacent tori. One set of parameters controls the waist radii, while the 

other set adjusts the heights of the centers, that is their signed distance from T. A degree 

of rigidity, in the form of matching conditions, is maintained to reposition the tori in response to 
the parameters, and the surface is smoothed using cutoff functions as needed. A single parameter 
^ works for [^, since there N = 2 and the symmetry between the sides of T forces ^ = 0. 

In the course of the construction it is necessary to solve for the proper selection of parameters 
along with the perturbing function. The dependence on the parameters of the “extended” com¬ 
ponents of the extended substitute kernel can be directly estimated with accuracy adequate for 
our purposes. It turns out that these components are primarily generated by dislocations resulting 
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from antisymmetric variation in pairs of ^ parameters associated to catenoids adjoining a common 
torus. 

The parameter dependence of the substitute kernel itself is more conveniently monitored indi¬ 
rectly, as in [^, via forces. The approximate kernel, and so the substitute kernel, on each torus, 
may be identified with approximate translations of the torus relative to T. In fact admits an 
exact Killing field which does not generate this variation of the torus but does approximate it lo¬ 
cally. The force in the direction of this Killing field through certain neighborhoods of a given torus 
then serves as an estimate of the projection of the mean curvature onto the approximate kernel 
and thereby as a proxy for the corresponding component of substitute kernel itself. The balancing 
equations and the analysis of the parameter dependence of the forces here are substantially more 
complicated than those of but no different in principle. 

Finally, estimates for the initial geometry, the linearized equation, the nonlinear terms, and the 
parameter dependence of the forces and dislocations are applied in conjunction with the Schauder 
fixed-point theorem to prove that given k, i, and N, for sufficiently large m there exists a choice 
of parameters and a smooth, appropriately symmetric perturbing function such that the resulting 
surface is minimal and a good approximation to the corresponding initial surface, so in particular 
embedded. 

Acknowledgments. This article presents the results of my PhD thesis, supervised by Nicos 
Kapouleas, who suggested the problems studied therein; I am deeply grateful for his guidance. 
Additional thanks are due to Scott Field, for discussions that further stimulated my interest in 
N > 2 stacking, and to Christina Danton, for help preparing the figures. 

2. Initial surfaces 

We realize as the unit sphere Z 2 ) G -t- |z 2 p = 1} in and set 

(2.1) T = |( 2 ;i,Z 2 ) : \zi\ = \z 2 

the Clifford torus whose equidistant axes are simply the coordinate unit circles Ci = {z2 = 0 } and 
C2 = {zi = 0 }. We define the covering map <!' : —)> §^\(C'i U C2) by 

(2.2) $(x, y, z) = sin (^z -h 0 cos (z + 0 ) , 

which maps (i) horizontal planes to constant-mean-curvature tori having axes Ci and C2, with 
<h({z = 0}) = T in particular, (ii) vertical lines to quarter great circles orthogonal to Ci, C 2 , and 
T, (hi) vertical planes of constant x to great hemispheres with equator C2, (iv) vertical planes of 
constant y to great hemispheres with equator Ci, and (v) vertical planes of constant x ± y to half 
Clifford tori through Ci and C2, orthogonally intersecting T along great circles. Writing gs for the 
round metric on and gE for the flat metric on we find 

(2.3) <^*gs = gE + (sin 2z) (dx^ - dy^) . 

The initial surfaces will be built by applying $ to a stack of horizontal planes connected by staggered 
catenoidal columns. 

Half-catenoids bent to planes. Since we shall make frequent use throughout the construction 
of cutoff functions, we fix now a smooth, nondecreasing T : M —?• [0,1] with T identically 0 on 
(-00,-1], identically 1 on [l,oo), and such that T — ^ is odd. We then define, for any a, 6 G M, 
the function 'ijj[a,b] : M —>■ [0,1] by 

(2.4) i/) [o, 6] = T o 
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Figure 1. An N = 5 initial surface, cut by four hemispheres and viewed in the 
(x, y, z) coordinate system; for clarity the picture grossly exaggerates the ratio, which 
in fact tends to 0 with large m, of the vertical spacing between the tori to the 
horizontal spacing between the catenoids. 


where La^b : M —)■ M is the linear function satisfying L{a) = —3 and L{b) = 3. Write B[{x, y),r, M^, 
for the open Euclidean disc in with radius r and center (x, y). Given z^^, zj’ G M and R, X,Y,t > 
0 with T < 2R < min(A, y), set 

(2.5) Tx,y,r = ([-y, y] X [-y, y]) \B[{0, O), r, m 2, ffij] 

and define the function (j)[zK, zt, R, X, Y, r] : Tx,y,t —^ Mby 


(I)[ZK, ZT, R, X, y, r] (x, y) =zk + (zy - zk)'iP [R, 2R] + y2^ 


( 2 . 6 ) 


+ sgn(zT - zk) arcosh ^ V' [2R, R] (^\/x2 + y2^ , 


where the sign function sgn : M —)• M takes the value 1 when its argument is nonnegative and the 
value —1 otherwise. Given further (xq, yo) G M^, define the embedding Text[{'^o, Yo); zjv, zy, R, X, Y, t] 
Tx,y,r M^ by 


(2.7) 


Text[{^o,Yo),ZK,ZT,R,X,Y,T]{x,y) = (xo + x,yo + y,(/)[zif,zr,i?,y,y,r](x,y)); 


and given also z)^ G M along with t' > 0, now assuming max(r, t') < AR < min(X, y), set 

( 2 . 8 ) 


Tx,y,.,y = ([-X,X] X [-y,y]) 


B 


X Y 


,9E 


UB 


X Y . -2 

) ,T ,9E 
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and define the embedding Tint[(xo, yo), R, X, Y, r, r'] : Tx,y,r,r' by 

(2.9) 


Tint{^,y) 


( {xo+x,yo + y,(j)[zK,ZT,R, f, |,r] (x + f ,y+ |^)) for (x,y) € [-^,0] x [-1^,0] 
< (xo + x,yo + y,(/) [z'^,zr,i?, f, |,r'] (x- f ,y - |)) for (x,y) G [0,X] x [0,y] 


[ (xq + X, yo + y, zt) elsewhere . 


The initial surfaces will be built from various applications of T^xt, for extreme or outermost 
tori and adjoining half-catenoids, and of Tint, for the intermediate tori and pairs of adjoining half- 
catenoids. The horizontal positions of the catenoids (equivalently the values of (xo,yo) in the 
parametrizations above) are decided by the lattice that the data k, £, and m set. The radii of the 
annuli of transition (determined by R) will be chosen on the order of but sufficiently smaller than 
the smaller of the two lattice edges. There remain N — 1 selections of waist radii (r and t') and 
— 1 more of waist heights {zk and Balancing conditions studied in the next section allow for 
the estimation of the values these 2N — 2 unknowns must assume for the construction to succeed, 
but their precise specification is made by the C and ^ parameters. Matching conditions then fix 
the heights (zt) of the tori, by requiring these to agree with the heights of the adjoining catenoids 
where they meet the tranisition annuli, in the case of the extreme tori, or to agree with the average 
of the heights of the upper and lower catenoids at the transition circles they adjoin, in the case of 
the intermediate tori. 


The hierarchy of data. For any positive integers k, i, and N > 2, the construction produces 
a sequence, indexed by m, of {2N — 2)-parameter families of initial surfaces. In order to obtain 
adequate estimates for (i) the initial mean curvature, (ii) the linearized operator, and (hi) the 
nonlinear terms, we will routinely make the assumption that m is sufficiently large in terms of k, 
i, N, and all of the parameters. Since we expect the ultimate parameter choices themselves to 
depend on m, it is necessary to assume that they are all bounded in absolute value by a constant 
c > 0 independent of m. Of course we do not yet know what range is needed for the parameters, 
but eventually we will be able to pick c in terms of k, £, and N so that for every sufficiently large 
m we will be able to find parameters bounded by c so that the corresponding initial surface can be 
perturbed to minimality. 

To continue with the definition of the intial surfaces we fix k, i, m, and N > 2 as well as c > 0 
and parameters C)C ^ \~c,c\^~^. For notational simplicity we assume k < £ and we write n for 
the greatest integer no greater than N/2, so that N = 2n when N is even and N = 2n + 1 when 
N is odd. We acknowledge a certain redundancy in the minimal surfaces ultimately exhibited, one 
which is easily removed by taking k and £ relatively prime. 

We associate to k, £, and m the rectangular lattice 


( 2 . 10 ) 



V2tt 

£m 


G 


i,j ez 


fixing the horizontal positions of the catenoidal tunnels. The tunnels connecting a given pair 
of tori have vertical axes either all intersecting this lattice or all intersecting its translate by 

( ^ ^ 0^ 

\ \/2km ’ \/2im ’ / ’ 

In the next section we will determine a collection of positive real numbers, as functions 

of k, £, and m, to serve as catenoidal waist radii when C = 0; the ratio of any two of these waists 
will be bounded independently of m, and we shall find 


1 


klmP" ( 1 
Ttt V 


£2 

2 


( 2 . 11 ) 


Tl = 


6 




Figure 2. A schematic profile of two N = A initial surfaces. The surface on the 
left has Cl = = 0, while the one on the right has nonzero values for all three 

of these parameters, with riCi = T 2 C 2 / U 3 C 3 , resulting in dislocation at the third 
torus from the bottom and nowhere else. (Of course the figure is not drawn to scale.) 


where C 2 is a positive number depending on k, i, m, and N but independently of m bounded away 
from 2. For general C we define the radii 


( 2 . 12 ) 


e'^^Ti for i = 1 

gCi+fc-T-im-2G^. for 1 < f < AT. 




(2.13) Zi = < 


n—1 ^ 

nCi - 2 ^ In - 2^t, In for f = 1 

j=i 1 

^ i—1 n—1 


+ 2 


i—lii — l+Tiii r)N mod 2^ 1 „ 1 

2 Z Tn m 

\j=l i=l; 

7V-1 

rjv-iCiv-i - 2^ In + 2 In for i = IV 

j=n 


In for 1 < i < iV 


heights of the catenoids’ centers are 

1 ( 


and the N — 1 
(2.14) 




Ti In 


Vi=i i=i. 


fi in- 

^ lOiniTj 


1 _ 2 N mod 2^^ - for 1 < z < iV ■ 

10£mTn 


These definitions can be understood as implementing the matching conditions mentioned earlier 
as well as the vertical offsets introduced by the C parameters. Each logarithmic term, ignoring any 
powers of 2 appearing as prefactors, represents the height achieved by a corresponding catenoid 
above its waist plane a distance from its axis, where the catenoids are meant to transition 
to planes (tori under 4>). The factor of 10^ is chosen—10 somewhat arbitrarily and i because we 
assume k < I —to ensure the transition is completed on the order (m“^) of the lattice spacing but 
well away from neighboring catenoids. Here the logarithm is used, as a matter of convenience, to 
approximate inverse hyperbolic cosine. 


7 















































Symmetries. We orient 
circle C in §^, we write 
>0 




by choosing the standard orientation on C^, and, given an oriented 


for rotation by 0 about C. Thus 


RC2^ZI,Z2) = 


{e^^zi,Z 2 ) and 


Z 2 ) = {zi,e''^Z 2 )- Note that <h intertwines each rotation of these two types with a hori¬ 
zontal translation. Let X denote reflection through the totally geodesic sphere with equator C 2 and 
(1,0) a pole, and let Y denote reflection through the totally geodesic sphere with equator Ci and 
(0,1) a pole. Then X(zi,Z 2 ) = (W, ^ 2 ) and y(zi,Z 2 ) = (zi,^); intertwines each of these with 
reflection through a vertical coordinate plane. 

We next define the subgroup Q of 0(4) generated as 


(2.15) 


G = 


2-7r 
Q km 

^C2 ’ 


27r 
^ im 


c. > X, Y 


Not only does Q preserve the lattice on T defined by k, i, and m, but it will also preserve the 
corresponding initial surface, whatever values the parameters assume, and we will later admit only 
functions equivariant under a natural action by G as candidate perturbations. Thus all elements of 
G will be symmetries of the minimal surfaces we construct. 

While intertwines reflections through a plane of constant x or constant y with an isometry of 
(reflection through the corresponding sphere), it is clear from 2.3 that (x, y,z) 1 —)> (x,y, —z) does 
not define a symmetry of but (x, y, z) 1 —)> (y, x, —z) does. Indeed reflection through a diagonal or 
antidiagonal line on {z = 0} corresponds to reflection through a great circle on T. Writing Z for 
reflection through the great circle {zi = Z 2 }, so that Z(zi,Z 2 ) = (^ 2 , 21 ), we see that intertwines 
Z with the above reflection through {x = y} on {z = 0} and further that Z preserves the lattice of 
the construction precisely when it is square. 

Accordingly, when k = i, we could enforce 


(2.16) 


Gs = 


277 
D km 

^Ci : 


X, ZR^f 


as the symmetry group of the construction. Actually, in the square case Gs will automatically 
preserve the resulting minimal surface provided we constrain the parameters to realize the sym¬ 
metries of Gs in the initial surfaces, cutting in half the number of free parameters. Specifically we 
would require that for N even Cn-i = Cn+i and ^n-i = —^n+i, while for N odd Cn-i = Cn-i-i-i-i and 
^n-i = —^n-i-i+ 1 ) leaving just N — 1 independent parameters. The entire construction, including 
the solution operator TZ, now coming with (N — l)-dimensional extended substitute kernel, will 
then respect Gs without any modifications. For this reason we will not again explicitly address the 
square case. 

Assembly and basic properties. Assuming the specification of the waist radii given in the next 


section and setting R = 
(2.17) 




\/2kn 


and Y = 


V2er 


let 


(Tea;i[(0,0),zf ,zi,ii’, A, y,ri](Tx,r,ri)) for i=l 
Gli = < ^ [Tint [(2f - 3) [f, X) ,zf_^,zf,Zi,R,X,Y,Ti-i,Ti] (Tx,y,r,_i,rJ) for 1 < i < A 
(Te,i[(2A - 2) (f , ,Z%_„ZN,R,X,Y,TN-l]iTx,Y,rN-i)) fo^ * = N. 

Then we define the initial surface 


N 


(2.18) 


Y = G\Jn,. 


i=l 


Proposition 2.19. Given positive integers k, I, and N >2, and a real number c > 0, there exists 
mo > 0 such that for every integer m > mo and every choice of parameters C)? £ [—c,c]^~^ the 
initial surface Y[N,k,i,m,C,f,], defined by 2.18, is a smooth, closed, orientable surface, embedded 








of genus kim?{N— 1) + 1, and invariant as a set under the action of Q[N defined 

by\KM 


Proof. In light of the selection of waist radii in the next section, especially 13.6| and 3.8, we have 

1 

lim Tj = lim In —--= 0, and so all of the claims are clear from the construction of S. □ 

m^oo m—>-oo lOrTTlTj 

Of course we have not yet precisely specified {Tj}, a gap filled in the next section. 


3. Forces and dislocations 


Forces. As mentioned in the introduction, we will eventually discover that on the toral regions the 
Jacobi operator C of the initial surfaces has approximate kernel, spanned by eigenfunctions with 
small eigenvalues. Actually we already have reason to anticipate this state of affairs. The Clifford 
torus T itself has Jacobi operator At + 4, with kernel spanned by cos \/2xcos \/2y, cos \/2xsin \/2y, 
sin- v/2n cos-v/2y) and sin \/2xsin-v/2y- Each of these functions is the normal projection of the 
restriction to T of a Killing field on S^. (The isometry group 0(4) of the sphere has dimension 
6, while its subgroup stabilizing the torus has dimension 2.) Evidently there are no exceptional 
Jacobi helds: all of them originate from families of ambient isometries. Note that while none of 
these four functions satishes the symmetries of our construction, the hrst one almost does, on a 
single fundamental domain containing (1,1) G T. 

The Killing field K inducing this Jacobi held is just rotation, toward C 2 , along the circle preserved 
by both X and Y. The circle has preimage under containing the z-axis, and with respect to the 
(x, y, z) coordinate system 

K = -^ cot fz + —sin V2x cos V^y H -= tan fz + —cos \/2x sin V^y dy 

( 3 . 1 ) \/2 V 4 ; ^/2 \ 4 :J ^ 

+ cos -v/2x cos \f2y ■ 


Thus K can be regarded as approximating clz, which of course generates vertical translations and 
corresponds to the constant 1 spanning the kernel of the Jacobi operator A for horizontal Euclidean 
planes—the inverse images under of constant-mean-curvature tori parallel to T. We will confront 
this basic eigenfunction, without any explicit mention of Killing helds, in the course of our detailed 
study of the linearized operator. Right now we intend to calculate the K force (also called flux in 
the literature) through various regions and to study its dependence on the initial surface. These 
forces will measure the projection of the surface’s mean curvature onto the approximate kernel, so 
in the hnal section we will apply the following computations to complete the proof of the main 
theorem. More immediately we will impose balancing conditions on the initial surface, such that 
the K force on various regions vanishes, at least within a margin on the order of the perturbations, 
by functions and parameters, that we will be making. This balancing will finally determine the 
waist radii, up to choice of C, thus completing the construction of the initial surfaces. 

Let T) = {K,r]i) , the K force exerted by flj, having outward conormal 77 *. One finds 


(3.2) 


T) -h 0{m ^Ti) 


+ l^zi for 7 = 1 

< 27r(ri - Ti-i) + for 2 < i < A - 1 

^ -2TrTN-i + for i = N, 


where the quantity 0{m~‘^Ti) absorbs terms which, for m sufficiently large—in terms of k, £, N, 
and c—are bounded by times a constant depending on just k, i, and N. Actually we are 

not in a position to justify this bound on the error terms given arbitrary {tj}. Our first application 
of the forces, however, is merely heuristic, as motivation for the definition of {t^} via 3.6 and 
3.4 for {ci} bounded independently of m. With this attitude one can simply retain the dominant 
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terms contributing to the force, captured on the right-hand side above, and impose the balancing 
described below to obtain |3.6| and |3.8[ Armed with these one can return to verify the estimate for 
the error terms encapsulated in the expression 0{m~‘^Ti). 

The computation of the forces themselves is simple. The term 27rrj corresponds to the force 
through a waist circle; the error originates from treating the waist as a Euclidean circle and can be 
assessed usingThe calculation of the force through a rectangular boundary component is eased 
by invoking the first-variation-of-area formula to replace the integral over that component by an 
integral over the exact constant-mean-curvature torus {z = zq} also bounded by it and having the 
same conormal there as the initial surface. This torus’ mean curvature is 2tan2zo5z5 accounting 
for the appearance of terms contributing to the force that contain height factors, and its area is 
27r^ cos2zo. The error here derives from approximating tanz by z and cos2zo by 1. 

Consequently Ti+i — Ti\- 0{m~‘^Ti) = 


27r(T2 - 2 ti) + ^ (2ri In for i = 1 and N>3 

27r(Ti+i - 2Ti + n-i) + ^ (in In for 2 < z < iV - 2 

27r(-2T^_i + TN- 2 ) + ^ {2tn-i In + giv-ir,v-i-giv- 2 riv- 2 ^ for i = iV - 1 and > 3, 


I'n + 0{m ^Ti) = 




Ti In 


1 


lOimri 


+ Tl 


27r(r„ - Tn-i) - 2^ (r„ In 


^0 


for re = 1 


1 


lOimTr. 


+ 




for re > 2, 


and Tn+i + 0{m ^ri) = 

2>r(T„+. - r„) + ((JV + 1) mod 2] ^ (r„ to + Mtot I- l T- n 

understanding = Cn+i = 0 in the last equation when N = 2. 


Balancing. The preceding equations will play an indispensable role in selecting the right param¬ 
eters at the conclusion of the construction, but at the moment we set C = ^ = 0 s-iid impose 
balancing so that all the above forces vanish up to terms of order rre“^r^, in order to determine the 
collection Starting with the last equation (possibly coupled with the one directly above 

it) and then proceeding inductively to compare the force through each with the force through 
its reflection (in the sense of Euclidean through {z = 0}, we evidently have 

. . Zn-i = Zn+i+i for N odd 

Zn-i = Zn+i for N even. 

Therefore, defining 


(3.4) 


Ci = for 1 < i < A" — 1, 


we need only solve the re equations 


27r(c2 - 2) = 


kim"^ 


2 In lO^rrer;^, provided re > 2 


(3.5) 


27r(cj+i - 2cj d-i) = 
27r{Cn-l-Cn) = 2^’^^^^ 


— —^2cj In lO^rezTj for 2 < i < re — 1 
ktm^ 


klm? 


Cn In lO^rrer^ (understanding cq = 0) 
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for {cj }”^2 'Ll- The first applicable equation directly yields 


(3.6) 


1 klrrfi /i C2\ 

r, = --- e 4,r I 2 > 

lOim 


where we understand C 2 = 0 when N = 2 and we note 3.3 implies C 2 = 1 when N = 3. By 
subtracting the other equations from multiples of the hrst we find 


C2C2 - C3 = 1 - 


Stt 


(3.7) 


kim? 

— Ci-l + C2Ci — Cj+l = — 


C 2 In C 2 for re > 3 
Stt 


klm? 


Cj In a for 3 < i < re — 1 


_ 2 (a^+i) mod - iV mod 2)cn = - 


Stt 

kim? 


, lnc„ 


Lemma 3.8. Given N > A there exist d 2 < < ■ ■ ■ < dn—depending on just N—with d 2 G (1, 2) 

strictly increasing in n, and there exists mo > 0 so that for each m > mo there are C 2 ,C 3 ,... ,Cn 
solving 15.71 and satisfying lim Ci = di. 

m^oo 


Proof. Bear in mind that balancing has been accomplished (by 3.6) for N = 2 and N = 3 (and 
in these cases the system 3.7 is vacuous). Momentarily ignoring the logarithmic terms, for N = A 
the system reduces to = 2, so ^2 = \/2 < 2, while for = 5 we get d^ — d 2 — A = Q, yielding 
^2 = < 2. Now the functions C 2 e-)• C 2 and C 2 e-)• — C 2 — 1 have nonzero derivatives at these 

respective values, so the the lemma is established for N = A and = 5 by applying the inverse 
function theorem and taking m large. Thus we may assume re > 3 and pursue an elaboration of 
the same strategy. 

For real f3 define the (re — 1) x (re — 1) matrices 


(3.9) 


and 


(3.10) 


A2niP) = 


i2n+l 


(/ 3 ) = 


//3 

-1 

0 

0 


0 \ 


-1 

/3 

-1 

0 


0 


0 

-1 

Id 

-1 


0 


0 

0 


-1 

Id 

-1 


VO 

0 


0 

-2 



[ (3 

-1 

0 

0 


0 

\ 

-1 

fd 

-1 

0 


0 

0 

-1 

Id 

-1 


0 

0 

0 


-1 

Id 

-1 

VO 

0 


0 

-1 

Id-lJ 


Thus, temporarily neglecting the logarithmic terms, we need to solve the system 


(3.11) An{^) 


(d2\ 

ds 


/1\ 

0 


subject to the constraints j3 = d 2 and d* > 0 for 2 < z < re. 



















where 


Using Cramer’s rule and expansion by minors 

Pn-i+i[N mod 2]{/3) 


(3.12) 


Pn[N mod 2]{P) 

Pi[0](^) = det^ 2 j(A) for i >3, 

Pi[l](A) = det A 2 j+i(A) for i > 3, 

P2[0](A) = A,P2[1] = A-1, 

Pi[0](A) = 2 , and Pi[l](A) = 1. 

Further expansion by minors reveals the recursive relations (independent of the parity of N) 

(3.13) P+i{X) = AP*(A) - P-i{X) for i > 2. 

On the other hand, the above expression for d 2 , applying the constraint d 2 = (3, can be rewritten 
as 

(3.14) PPM=Pn-im, 
whence 13.131 delivers 

(3.15) P„+i(/3) = 0. 

We now claim that for each n >3 (and either parity of N) 

• Pn has a root strictly greater than 1 ; if 7 n is its largest such root, then 

• Pn-i{x) > 0 whenever x > 7 ^, and 

• Pn+l{'yn') ^ 0 . 

Moreover, 7 ^ is strictly increasing in n. 

These claims can be established by induction on n. The case n = 3 is easily verified: ^ 2 ( 0 ] (a:) = x, 
^ 3 [ 0 ](x) = - 2, and P4[0]{x) = x^ - 3x, so 73 ) 0 ] = y/2, P 2 [ 0 ](x > 73 ( 0 ]) > 0, and P 4 [ 0 ]( 73 [ 0 ]) = 

2^/2 — 3y/2 < 0, while P 2 [l]{x) = x — 1, P 3 [l](x) = — x — 1, and P 4 [l](x) = x^ — x^ — 2x + 1, 

so 73 ( 1 ] = P 2 [l](a: > 73 [1]) > ^ 2 ~^ > 0, and P 4 [l]( 73 [l]) = < 0. Now suppose the 

claims hold for n = j. According to the third claim Pj_|_i( 7 j) < 0, but Pj+i(x) > 0 for large x, so 
this implies that Pj+i has a root greater than 7 ^, so_ 7 ^i exists and 7^+1 > 7 j. Therefore Pj(x > 
7 j+i) > 0 by the maximality of 7 j. Finally, using |3.13| Pj+2(7i+i) = 

which is negative, since the first term vanishes and the second has just been established positive. 

Thus /3 = 7 n+i solves 3.15 Since 7 ^ is greater than 1 and strictly increasing in n, we find, using 
again [3T3 that 

(3.16) Pjipin) — Pj+l(7n) = (f ~ ln)Pj{'^n) + Pj-l{ln) > 0 

whenever j < n. It then follows from |3.1^ that d 2 < d^ < d^ < ■ ■ ■ < dn- 

Finally we claim d 2 <2. In fact we assert that for each n >2 (regardless of the parity of N) 

• Pn{x) — Pn-i{x) > 0 and Pn{x) > 0 whenever x > 2 , 

which is proven by induction on n. For n = 2 and x > 2 we have P 2 (x) — Pi(x) = x — 2 > 0 
(whatever the parity of N) and clearly both P 2 [ 0 ](x) = x > 0 and P 2 [l](x) = x — 1 > 0. Assuming 
then that the claim holds for n = j, we get from ]3.1^ assuming still x > 2 , 

P,+i(x) - Pj{x) = xPj{x) - Pj-iix) - Pj{x) = (x - l)Pj(x) - Pj-i(x) 

> Pj{x) - Pj-i{x) > 0 


(3.17) 


and therefore Pj{x) > 0 as well. We conclude that for every n > 2 we have Pn{x) > 0 whenever 
n > 2, so all roots of Pn lie to the left of 2, establishing the bound on d 2 - 
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It remains to reintroduce the logarithmic terms. We define the function Fn 
(3.18) 


pn—1 


pn—1 


by 



( X2fi 







xz 



3^3 



Fn 


= An{x2) 





\XnJ 



\XnJ 



m=2- 










fd2 0 0 

• •• o\ 





dz 0 0 

... 0 


— 

^i)i = 2 

^Af(<^2) + 








, 

0 ■ 

0 ■ 

••• V 


( d2 

-1 

0 

0 



0 \ 

dz 

13 

-1 

0 



0 

d^ 

-1 

(3 

-1 



0 

dn—l 

0 


-1 

13 

-1 

y dn 

0 


0 - 

2(A^+1) mod 2 

(3 — N mod 2y 


(3.19) 


whose determinant is 


det A]sf{d 2 ) + det 

(3.20) 


=Pn(d2) + X: i+1 (^ 2 ) 0. 

i=2 

We conclude by invoking the inverse function theorem and taking m large. 


□ 


Parameter dependence. Henceforth we fix {ci}”^2 in the previous subsection, the corre¬ 


sponding Zi as defined by 3.6 = CiZi for 2 < i < n, and for remaining values of the index in 


accordance with 3.3 The next lemma completes our estimation, at this stage, of the forces. Both 
the C and ^ parameters influence the forces, which will be analyzed to address the substitute kernel. 
The ^ parameters also control dislocations {Vi}, or displacements of a pair of adjacent inequivalent 
catenoidal regions relative to the toral region they share. 

Lemma 3.21. Given c > 0 an d positive integers and N >2, parameters Ci? £ , 

and {Ti}f^^ as defined by 2.12. let 


(3.22) 


Vi = 


_ I \Tiii - \Ti-iii-i for2<i<N -I 


0 for i = 1 and i = N. 

Then for m sufficiently large in terms of k, i, N, and c, 

(i) ^{J^i - Vi+i) + ^{Vi + A+i) + 0(1) = {ZQi, where 
(a) the {N — 1) X (N — 1) matrix 


(3.23) 


Z = 


/ 87r 

-C2 

0 

0 


0 


87rc2 

1 + C3 

-C3 

0 


0 

0 

8 vrc3 

-C2 

C2 + C4 

-C4 

0 


0 

87rc4 

0 

-C3 

C3 + C5 

-C5 


0 

87rCN-2 

0 


0 

—CN-Z 

CN-Z + CN-1 

— CN-l 


\87rcAr_1 0 


—CN-2 CN-2 J 


0 
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(understanding Cj = 0 for j > N and j <0) has inverse bounded independently of m and c; 
(Hi) + Tn) = Cl + CN-if.N-1 + 0(1); 

(iv) ^Vi = Cifi - Q-iCi-i + 0(1); and 

(v) the {N — 1) X {N — 1) matrix 


(3.24) 


/ 1 0 0 ••• 0 

— 1 C2 0 0 

0 — C2 C3 0 

0 0 — C3 C4 0 


0 cn-1^ 

0 0 

0 0 

0 


0 0 
\ 0 0 


0 —CN-3 CN-2 0 I 

0 0 —CN-2 CN-lJ 


has inverse bounded independently of m and c. 


Proof. The calculations are straightforward, using our previous relations between the forces and 
taking advantage of the balancing just imposed for () = C = 0. As for the invertibility of the 
matrices, expansion by minors and a simple inductive argument reveal that the determinants are 
strictly positive provided merely that each Cj is strictly positive, which condition is obviously 
satisfied; that the inverses are bounded as claimed follows from the convergence of c* to di with 
large m. □ 


4. Estimates of the initial geometry 


Norms and spaces of sections. To state the estimates and to carry out the rest of the con¬ 
struction, we must first identify certain norms and corresponding spaces of sections. For the most 
part our notation is standard and speaks for itself. Given a vector bundle E —)• M, a nonnegative 
integer j, and an exponent a G (0,1), we write and T\'^^{E) for the space of sections of 

E having component functions of class or (7/^“ respectively relative to every local chart and 
trivialization. We set T°^{E) = H^o ™ ^ trivial bundle M x M, we denote 

the spaces here and below by replacing T with C and E with M, as usual. 

All of the vector bundles of interest to us are derived from tangent bundles by a combination 
of duality, tensor product, pullback, and projection; a single Riemannian metric will determine 
canonical metrics and connections on all these bundles. When there is no danger of confusion we 
write simply H-H for the corresponding pointwise norm and D for the connection. Given a section 
u of a bundle E —)■ M so equipped, we define the standard global norms 


j 

(4.1) ||u ; W(E,5r)|| = ||u||. = V sup \\D^u{p)\ 

i=o 


as well as the Holder seminorms 
(4.2) [u]a = sup 

7:[0,1]-5>M 


«( 7 ( 1 )) - ^oM7]«(7(0)) 

i7r 


where the supremum is taken over all piecewise paths and, given such a path 7 ; [a, 6] —)• M 
and parameters s,t £ [a, 6], ^*[7] : —)■ is parallel transport along 7 from the fiber over 

7(5) to the fiber over 7(f) and I7I denotes the length of 7. 

Then we can define also the Holder norms 



Note that for functions on convex open subsets of Euclidean space these Holder norms agree with the 
conventional ones, and on precompact subsets with sufficiently regular boundary the two norms are 
at least equivalent. Generally, the spaces and T^{E,g) consisting of sections with finite 

corresponding norm enjoy many of the properties familiar from the Euclidean case; in particular 
T^'l^{E,g) embeds compactly in T^'^{E,g) whenever M is closed and 0 < a < /? < 1. 

In the construction we will routinely wish to compare norms of the above type induced by 
different metrics on a single manifold. The definitions make it easy to see that ||n : r^’°‘{E,h)\\ < 
C\\u : r^'^{E, g)\\, where C is controlled by the gf-norms of h, its inverse, and finitely many g- 
derivatives of h (the maximum order needed depending in a transparent way on j, a, and on the 
bundle E). 

If M is a two-sided hypersurface immersed in a Riemannian manifold N, and Qm is a group 
of isometries of N preserving M as a set, then Qm acts on a section u of the normal bundle of 
M by {q.u){p) = (0~^(p))] for each g G Qm] because this bundle is just the trivial M bundle 

over M, its sections can be identified with functions (for us representing mean curvature or normal 
perturbations) on which the corresponding action of Qm is given by (g./)(p) = (—1)®/ (g”^(p)), 
where (—1)® is 1 if g preserves each of the two sides of M and —1 if it exchanges them. 

All the elements of the symmetry group Q of the construction can be seen to fix each side of 
S, so this action is trivial, but we remark that Qs contains reflections through great circles which 
when N is odd exchange the sides. As already explained though, our construction will produce 
^g-invariant minimal surfaces from ^^-invariant initial surfaces despite enforcing just Q throughout. 

In general we will append the subscript Qm to a space of functions to denote that subspace 
consisting of function which are equivariant under the Qm action just described. Sometimes we will 
append the subscript c to a space of sections to isolate those sections having compact support. 

Finally we will often wish to work with weighted versions above the norms. For this construction 
the following definition suffices; 


(4.4) 


u:r^’'^{E,g,f)\ 


sup 

peM 


Up) 


where f : M ^ (0,oo) is a weight function and B{p, l,g) is the ball of sf-radius 1 centered at p. 
We will also make use of weighted E-^ norms, with the obvious definition. 


The y metric. It is the primary task of this section to estimate the intrinsic and extrinsic geometry 
of the initial surfaces. To fix the extrinsic quantities we pick on each initial surface S the global 
unit normal u which is directed toward Ci at the points of S closest to Ci. We then define A to be 
the scalar-valued second fundamental form of S relative to i', and we take the scalar-valued mean 
curvature H to be the trace of the former. 

Every initial surface admits by virtue of its construction a natural decomposition into overlapping 
regions, each of which resembles either a portion of a torus or a truncated catenoid. Modulo the 
horizontal symmetries, there are N such toral regions, one for each torus incorporated in the 
construction, and there are — 1 catenoidal regions, one for each pair of adjacent tori. Definitions 
are made in the subsections below. The estimates will then be obtained by treating the catenoidal 
regions as perturbations of Euclidean catenoids and the toral regions as graphs over the Clifford 
torus. 

Because all these regions shrink with increasing m and because even on a fixed initial surface the 
characteristic scale m~^ of the toral regions dwarfs the characteristic scale ri near the waists, it 
will be advantageous to uniformize the problem by working with a metric x on each initial surface 
conformal to the natural one g = i*gs induced by the round spherical metric gs and by the inclusion 
i of the initial surface in We will set 

(4.5) X = p^g, 
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where the conformal factor p~^ is to be defined as a ^-equivariant function on S (i) measuring 
on each catenoidal region the ^*~^gE distance to the axis and (ii) transitioning smoothly to the 
constant m~^ by the edge of the toral regions. 

To be precise, if we first define /o[zi, Z 2 ] : —)> M, for given zi, Z 2 G M, by 


(4.6) 


p[zi,Z 2 ](x,y,z) 



0 otherwise 


_^1 

lQlm\ 




if Zl < Z < Z2 


and p : —)• M by 

(4.7) 

N-l 

p(x, y, z) = m+ ^ zf+i] 

(p>9)ez2 


{2p + (i + 1) mod 2)7r 
X jY' 


(2g + (i + 1) mod 2)7r 


V2k 


m 




m 


taking z^ = — | and = |, then p G (7^(2) is uniquely defined by 
(4-8) = ^$-i(s). 

Equipped with the x metric, each catenoidal region tends with large m to the flat cylinder of 
radius 1, while each toral region tends, away from the catenoids adjoining it, to a flat x 
rectangle. 

Before proceeding, we briefly mention a couple differences of our approach from |^. First, our 
catenoidal and toral regions above correspond to their extended standard regions, but since we 
never view their standard regions or transition regions in isolation, we omit the modifier extended. 
Second, our use of the x nietric follows theirs to study the mean curvature equation on the initial 
surfaces globally, but whereas Kapouleas and Yang introduce another metric, h, conformal to g in 
order to analyze the approximate kernel, we will apply the x rnetric to this problem as well, in the 
next section. 


Catenoidal regions. Given a > 0 set = [—a, a] x and, given further (xo,yo,zo) G and 
T > 0, define the embedding k[(xo, yo, zq), t, a] : by 

(4.9) k[(xo, yo, Zq), t, a]{t, 6) = (xq, yo, zo) + r(coshf cos 6, cosh t sin0, t). 

If we endow the cylinder M x with its flat metric 

(4.10) xk = dt^ + 


and the image of k above with its natural catenoidal metric g induced by the ambient Euclidean 
metric gE on the target, we can easily calculate that k is a conformal map with 

(4.11) k*g = (r^ cosh^ t) xk- 

We also observe that the pointwise inner product of the Killing field on with a global unit 
normal for the catenoid defines a Jacobi field. Multiplying the Jacobi equation for this field by the 
conformal factor sech^ t one can readily solve for the squared norm 


(4.12) 

of the second fundamental form of k. 
Next, for 1 < i — 1 set 



2r ^ sech"^ t 


(4.13) 


Oj = arcosh 


1 

lO^mr* 
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and define k* : Ka^ —>■ by 


Then the image of Ki is entirely contained in the initial surface and defines the associated catenoidal 
region 

(4.15) )C[i] = Ki . 

Since the covering map is, on small scales, an approximate isometry, we expect that each 
catenoidal region will, in a suitably rescaled sense, converge in the large-m limit to an exact 
catenoid in Euclidean space. The next proposition quantifies this convergence. 

Proposition 4.16. Given positive integers k, I, and N >2, along with c > 0, there exists mo > 0— 
depending on k, £, N, and c—and there exists C > 0 —depending on k, I, and N, but not on c—such 
that for each [—c, c]^~^, for each m > mo, and for 1 < i < N — 1 

(^) 

(ii) 

(Hi) 


-Xk-.T^ xk) II < Cm^n; 


ty* ^—2 i 


9S 


— 2sech^t : (]Ka.,y;x) < Cm and 


Ip {Ka,,XK,'inn?K*p +m?Tl)\\ < C. 


Proof. Clearly K*p = ^secht, and from 


4.14 


4.9 


and 2.3 it is simple to calculate 


(4.17) 


CiX~ Xk = sin 2z (tanh^ t cos 20 dt^ — 2 tanh t sin 26 dt dO — cos 26 dO^) 


where z = + Tit; since xk = dC + dO^, since dz = Ti dt, and since the definitions of r^, zh, and 

Oi imply |z| < Cm^ri, this establishes (i). 


It is fairly efficient to estimate 


by the same method outlined above to find 


A 


. In this 


context too, that is relative to ^*gs, dy_ is a Killing field, so / = {v, d^) defines a Jacobi field on Ki. 
We find 


(4.18) 

with normalization factor 

(4.19) N = 


N sech t cos 0 „ N sech t sin 0 „ 

V = -:-—(9„- ^—dy + N tanht 


1 + sin 2z 


1 — sin 2z 

cos 2z 


\/1 — sech^ t cos 20 sin 2z — tanh^ t sin^ 2z 


f{t, 6) = N sech t cos 0. 
is just 2gs, we have Jacobi equation 


and so 

(4.20) 

Since the Ricci curvature of 

(4.21) 

on ]C[i]. Multiplying by p~‘^ and rearranging, we obtain 

d!f + dlf 1 


\f + 


/ + 2 / = 0 


(4.22) 


-2 


P 


f 


j (^X ^x) f ‘^p ’ 


which delivers (ii) in light of 4.17, the expressions just given for / and N, the estimates noted above 
for z, the inequality p > Cm, and the equality dp~^ = ptanhtdt on }C[i]. 

Finally, by (a) considering the family of immersions {^^[s] = Ki + where + incidates the 

obvious application of the Euclidean exponential map transferred via <I> to a neighborhood of S 
in (b) calculating the rate of change ^|^_q at s = 0 of the volume form of the family 

17 
























in {t, 9) coordinates, and (c) appealing to the first-variation-of-area formula, one can without too 

much labor compute the exact mean curvature 

(4.23) 

tanh^ t sin 4z cos 2z + ^Tp‘^[l + tanh^ t] cos 29 sin4z + sech^ t tanht[4 sin 2z — (1 + 3 sin^ 2z) cos 29] 

(1 — sech^ t cos 29 sin 2z — tanh^ t sin^ 2z)^/2 

proving (iii). 

□ 

Graphs over immersions. The estimates away from the catenoidal regions will be obtained by 
treating the initial surface there as a graph over the torus, as an application of the following lemma. 

Its last two items will be used again in the final section to estimate the contribution to the mean 
curvataure of the perturbed surface which are nonlinear in the perturbing function. 

Lemma 4.24. Let (j) : T, ^ M be a smooth two-sided immersion of a surfaee S into a eomplete 
Riemannian 3-manifold M with smooth metric g, and fix a global unit normal v G cjfTM for f. 
Given u G define the map : S —)■ M by fuip) = exp^^^^ u{p)v{p), where exp : TM —)■ M 

is the exponential map on {M,g), and define the function ; S —)■ M 

('‘■25) = E (l|o''.4(f')IF + l|o'fiUfe|.*)|,,,||r) . 

j=0 ^ ^ 

where B{p, \u{p)\ ,g) is the open ball in M of g-radius \u{p)\ centered at p, A is the second funda¬ 
mental form off relative to u, R is the Riemannian curvature of{M,g), and each instance of D 
is the appropriate connection canonically induced by g and f. 

Given a G (0,1), there exist constants C,e > 0 —with C depending on a but both independent of 
M, g, S, f, and u—such that if \\uJuU : C'‘^(S)|| + \\du : T® (r*S, </>*g')|| < e, then fu is a two-sided 
immersion satisfying the following estimates, where covariant derivatives, norms, and seminorms 
are defined relative to ({Ag, where Au and Hu denote the second fundamental form and mean cur¬ 
vature of fu relative to the global unit normal Vu distinguished by the condition that Vuip) have 
positive inner product with exp^^^^ and where the Jacobi operator C is defined by 

^ sL=o 

(i) \Wu9 - f*g\\ < C \u\ + 

(a) WDflgW < C {ujI |ii| -\-LVu \\du\\ + ||du|| ||L>2 u||); 

(iii) [Dflg]^ < C (^||w„||o [du]a + ||du||o + ||a;„||^+" ||tt||o 

+ ll‘^^^llo Iloilo + \\du\\o + [du\a Hd^wIIq); 

(iv) \\Au - A\\<C (ll-D^nll +oju \\du\\ -\-ujf |tt|); 

(v) [Au - A]^<C ||u||g + ||a;„||o+“ ||(iu||g + ||a;„||o [du]a + IM^^IIo [du]a 

+ ll^^^llo \\du\\l ll-^^^llo + 

(vi) \\Hu — Hq — Cu\\ < C |ii|^ + uju |it| ||-^^^|| + 1 

(vii) [Hu - Hq- Cu]^ < C ||u||q + ||a;„||^+“ ||tt||o ||du||Q + ||a;„||o ||ii||o [du]a 

+ lk«llo’^“ \\du[[l + llw^llg ||(iii||o [du]a + lldullg [du]a ll-D^^^llo 

+ ||w„||o ||dii||o ||L>2u||g + ||dn||o [D‘^u]^ + ||cj„||o [u]a ||-D^«||o 

+ Ikwllo’^^ ll'“llo l|■^^'“llo + Iloilo \d'<Aoi ll-C'^'^llo 

+ Ikwllo ll'^llo Ma + ll‘(^«llo Iloilo 
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We omit the proof, which consists of elementary calcnlations, making repeated use of the Jacobi 
equation along geodesics normal to S. 

Toral regions. The overlap of the toral and catenoidal regions is limited by a collection 
of positive constants, bounding below the absolute value of the t parameter on the corresponding 
catenoidal regions and chosen large a bit later in ter ms o f k, i, an d N, but bounded independently of 
m (unlike each at, which is of order m?). Recalling 

and R = define Ti ; Tx,y,ri cosh6i by 

(4.26) Ri(x,y) = $ {Text [(0,0), zf, zi, R, X, T, n] (x,y)) , 

define Tjq : Tx,Y,rjv_i cosh6jv-i by 

(4.27) Tn{^, y) = $ {Text [(iV - 1)(X, R, X, Y, r^_i] (x, y)) 


2.7 


and 


2.9 


and keeping X = 


Y2kv 


Y = 




and for 2, Y i ^ N 1 define 7) . Tx,y,ri_icoshfei_i,Ticoshfti 
(4.28) Ti{x,y) = ^ (Xint 


§3 by 


(2f - 3) ( y, ^ ) ,zf_^,zf,Zi,R,X,Y,Ti-i,Ti 


(x,y) 


The image of R, entirely contained in the initial surface, gives the toral region 

(4.29) r[f]=imri. 

The parametrizing rectangles naturally carry the flat metric gE = dx^ + dy^, but we equip each 
one also with the conformal metric 

(4.30) XT = {Ttp^) qe. 

In the next section we will define the extended substitute kernel needed to complete the con¬ 
struction, as outlined in the introduction. Then, in the final section, the role of the dislocations 
will become clear; the dislocation Pj (recall 3.22) on the toral region T[i] will be varied to cancel 
the “extended” portion of the extended substitute kernel supported there. For this reason it is 
necessary to isolate the dominant contribution of each dislocation to the mean curvature, and to 
that end, for 2 < i < X — 2, we define Vi G C°°{T[i]) by 


T*Vi{x,y) =•!/> 


1 1 

5£m’ lOim 


(4.31) 


- V’ 



1 


5£m 10£m 


and Wi G C^{T[i]) by 
(4.32) 


Wi = p 


rp-l* Vi 

T 9 e * 


Restricted to T[i], the function Vi should be regarded as the section of the normal bundle gen¬ 
erating dislocations there, and we will see the function Wi then captures the principal effect of 
dislocation on the mean curvature; later the collection {wi}^^^ will reappear as the defining basis 
for the extended part of the extended substitute kernel. For notational convenience, we define 
wi G C°°(T[1]) and we £ C°°{T[N]) to be identically 0. 

Proposition 4.33. Given positive integers a, k, I, and N > 2, as well as c > 0, there exists 
mo > 0 —depending on k, £, N, and c—and there exists a constant C > 0 —depending on k, £, and 
N, but not on c—such that for each C, G [—c, c]^~^, for each m > mo, and for 1 <i < N 
(i) \\Tfx - XT : F^’" (T*7;-V[i]®2, xr) || < (mVi + sech^ min^- bj); 
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(ii) 

(in) 


'* ^—2 


: (7°’" (T^ ^7~[i], xt) < C [m ^ + sech^ miiij bj); and 


Ttp 

p~‘^H — ViWi : C°’“ (T[i], X) rn?p~‘^Ti + m^rf) || < C*. 


Proof. First observe that if T[i] n IC\j] / 0, then 


(4.34) 


* — 1 * 

XT 


T[i]nK[j] 


= tanh^tdt^ + dO"^, 


which with item (i) of 4.16 proves item (i) of the present proposition where the toral and catenoidal 
regions overlap. Similarly (ii) and (iii) have already been established on these intersections. 


We will finish the proof by verifying the estimates on {p < 100m}. For this we apply 4.24 
viewing Tj as a perturbation of the embedding cj) = w o Ti oi inside the Clifford torus 

T, with t*7 : —)• T indicating orthogonal projection in onto T. As the function generating the 
perturbation (that is playing the role of u in the statement of 4.24) we have, for 2 < z < — 1, 

Ui ■ TjfjyjTi-i coshcoshbi ^ given by 

1 1 


(4.35) 

where r±i = 

(4.36) 


Uj(x, y) =Zj + ( zf_^ + Ti-i arcosh ^-z, 

'^i—1 


if 


+ 


zf — Ti arcosh — — Zj 1 -0 


1 


5£m’ lOfm 


5£m’ 10£m 

^ (n), 


(r-i) 


\/(^^ + ^ M given by 


2V2k 

Ul = Zi + 


zf — Ti arcosh — — zi 

n 




Mm’’ IQim 

where ro(x,y) = s/y? + y^; and, for i = N,un ■■ Tx,y,rjv_i coshfejv-i ® 


(ro), 

given by 


(4.37) 


UN = ZN + [ Zn -1 + tn -1 arcosh 


rp 

TN-1 


-ZNjl) 


1 


1 


5im Wim 


(rp). 


For 1 < i < A^ and m sufficiently large in terms of k, £, N, and c, we estimate 

(4.38) \\ui : {T~^ iT[i]\{p < 100£m}) ,m‘^gE)\\ < Cm^Ti, 

where we emphasize that C is independent of c and m. Thus |4.24| yields 

(4.39) \\n?T*gs - m^gE : F^’" (T*Tr^ (r[f]\{p < 100£m}), m^^g) || < CmVi. 

T* o 

Next observe that and its reciprocal are uniformly bounded on {p < 100£m}, independently 
of m (and c), as is the mfgE norm of every wfgE covariant derivative of Thus every Holder 
norm defined by xt is controlled, uniformly in m and c, by the corresponding norm defined by 
mfgE, securing (i). Similarly, since T has s econd fundamental form (dy^ — d'x^)dz relative to the 
coordinates induced by by appealing to 4.24 and comparing xr to mfgE on {p < lOOmj, we 
obtain (ii). 

Finally, the Jacobi operator on T is just + 4, and we find 


A 


qe'^i 


— Vip^Wi : (7°’“ (r^ ^T[i] n {/9 < lOO^mj, vnf'gE) || < Cm^Ti. 


(4.40) 

(Note that without subtracting the dislocation term, this contribution from the Laplacian of Ui 
to the mean curvature could not be bounded independently of the parameters.) With this last 
estimate now in place, the proof is completed by another application of 4.24, the same comparison 
of XT norms to mfgE norms, and the equivalence (by constants independent of m and c) on every 
subset of T[f] of the (7°’" norms induced by xr and Xi which itself follows from (i). 

□ 


2P 




































Decay norms and a global estimate of the mean curvature. As mentioned in the intro¬ 
duction, because the characteristic size ri of the catenoidal waists is so much smaller than the 
characteristic size m~^ of the toral regions, we must allow perturbing functions to be much larger 
on the toral regions than on the core of the catenoidal regions. For this reason we will weight 
our norms by powers of the factor mp~^, which takes the value 1 a maximal distance from the 
catenoidal regions and has order mri at the waists. Specifically, for each a G (0,1), 7 G [0,oo), 
and nonnegative integer j, we dehne 


(4.41) 

and set 

(4.42) 




: C^’' 


m ' 


S,X, 






7>a.7 


< 00 


}. 


where, as explained in the first subsection, the subscript Q has the effect of admitting only functions 
invariant under Q. 

In order to secure acceptable decay estimates for solutions to the linearized equation, we will 
need the following estimate for the initial mean curvature. 


Corollary 4.43. Given a G (0, 1) and positive integers k, i, and N > 2, as well as a real number 
c > 0, there exists a constant C > 0 —depending on k, i, and N but not on c—and there exists 
a positive integer mo—depending on k, I, N, and c—such that for each m > mo, each ^ 
[— c, and each 7 G (0,1) 

N-l 

^=2 0 , a ,7 

Proof. The estimate follows from the definition 2.12 of ti, from the definition |4.8| of p, from item 
(iii) of 4.16[ using also the equivalence of y and xk norms implied by item (i) of 4.16 and from 


(4.44) 


iW, 


< Cti. 


item (iii) of 4.33 


□ 


5. The linearized operator 


To review the strategy outlined in the introduction, given an initial surface S, embedded in 
by i : S —)■ S^, and given a function u G we define the perturbation : S —)■ by 

^'u{p) = exp^(p) u{p)v{p). As observed in 4.24 when u is sufficiently small, iu will be an immersion 
with well-defined scalar mean curvature relative to the unit normal obtained as a perturbation of 
p. For fixed k, I, and m, the mean curvature of may be considered a functional T-Lc^^^lu] of u 
as well as of the parameters f and f. We intend to solve the equation = 0, when m is 

sufficiently large in terms of k and I, and to obtain estimates for u demonstrating that Lu is a 
smooth embedding close in a precise sense to the original l. 

A major step toward the solution consists in the study of the initial surface’s Jacobi operator, 
the linearization C ofPL defined by 


(5.1) 


Cu= ^ 
dt 


t=o 


PL[tu] = (^Ag + 


-\- 2] u 


Actually, because of the uniformity afforded by the y metric, it is much more convenient to study 
instead 


(5.2) 


^x~ P P 




from 


4.16 


and 


4.33 


which, recalling the estimates of p ^ 

(S) —)• {Td) which is bounded independently of m and c. 


/^2,a,7 


clearly dehnes a linear map 
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In this section we construct a likewise bounded right inverse 7^ to modulo extended substitute 
kernel. We do this by first analyzing “semilocally”, meaning on the toral and catenoidal regions 
individually, and by observing that on each of these has a simple limit as m —)■ oo. We 
invert these regional limits (modulo extended substitute kernel in the toral cases) and so produce 
approximate semilocal inverses to which are applied iteratively, using decay properties of the 
solutions they yield, to construct TZ. 


Continuity in the parameters. Because we will have to select the parameter values of the initial 
surface in parallel with the perturbing function, it will be necessary to compare functions defined on 
initial surfaces corresponding to different parameter values. Such comparison will be accomplished 
by pulling back functions through certain diffeomorphisms between the initial surfaces. We define 
these diffeomorphisms as compromises between natural identifications on the various standard 
regions. More precisely we set 

(5.3) a,- = arcosh —-- 

lOimZi 

so that Oj is simply the value taken by ai when C = 0. Then we find that the functionTj : [r^, oo) —)• 
ti, oo) given by 


(5.4) rj(r) = Zi cosh ( — arcosh — ) ■0 

d'i TV 


1 


1 


20im ’ 30£m 


(r) + r^l; 


1 


1 


30im’ 20im 


(r) 


is the identity on oo) and strictly monotonic everywhere for m sufficiently large in terms of 

k, i, N, and c. Now we can define the diffeomorphism 


(5.5) : S[A1, k, £, m, C, t ^[N, k, i, m, 0, 0] 

by demanding that it commute with the action of Q, that 

^ o ,^](t, 9) = Hi[£), 0] ((sgnt) arcoshr~^r^j(ri coshf), 0) , and that 
^C.?lr[i]\(K:[*-i]uK;[i]) ^ o Ti[C,^] V 

understanding /C[0] = IC[N] = 0. Obviously each diffeomorphism preserves the catenoidal 
regions, but we arrange also for preservation of the toral regions by decreeing 

(5.7) h = -b 

Oii 

with h a positive constant at our disposal but independent of m and the parameters. Finally we 
set 


( 5 - 8 ) = Pip 

Lemma 5.9. The map 
(5.10) 

from the set of parameters [— c, x [—c, to the space of bounded linear maps from (S[C = ^ = 0 ]) 

to (7°’“’''' (S[C = = 0]), equipping this last space of maps with the operator norm, is continuous 

with bounded image. 

Proof. In fact it is obvious from the definitions of Cy. and that the coefficients, relative to 
either {t,6) or (x, y) coordinates on ^ = 0], of the operator in question are smooth in the 

parameters and coordinates and are uniformly bounded in (^^(^[(^ = ^ = 0],y). □ 
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Catenoidal solutions. From items (i) and (ii) of 4.16 we can see 

c-l 


(5.11) 

We set 

(5.12) 
and 

(5.13) 

and define : C, 


lim K*CyK* 

m^oo * * 


= Ck ■= ^XK + ^ sech"^ t. 


2,q,7 


m) 


(5.14) 


Gk.[{\ = {0 e ^ I g^i] c /C[f]} 

C^’"’^(/C[z]) = C'^’“J/C[z],x,mV 

C°’“’"(/C[z]) by 


-7'! 


and in the next proposition construct a suitable inverse. 

Proposition 5.15. Given 0,7 G (0,1), c > 0, and positive integers k, i, and N, there exists 
mo > 0 —depending on the foregoing data—and there exists C > 0 —independent of c—such that 
for C, ^ G [—c, I < i < N — 1, and m > mo, the map 

(5.16) ■■ (^Wc=?=o) ^ (/C[z]c=s=o) 

has norm bounded by C, is constant in and depends continuously (with respect to the operator 
norm) on (, and moreover there exists a linear map 






m) 


(5.17) Uk 

such that if f & {^[^]) = ^A:[i]/) then 

(i) ^K[i\u = f; 

(ii) \M 2 ,a,x < <^ll/llo,a, 7 / ^^^d 

(Hi) the map {lC[i]c_=^=o) 6 '^’“’"’'(/C[i]^=g=o) has operator norm bounded 

uniformly in the parameters and m by C, is constant in f,, and depends continuously (with 
respect to the operator norm) on (. 

Proof. Set K = M X and define 

(5.18) C^’"’^(K) = juG | u{t,e) =u{t,-e) = n(t, 7 r- 0 )|. 

By separation of variables we can reduce the study of Ck '■ C‘g°^''^(K) —)• C'g’“’"’'(]K) to the study of 
the ordinary differential operator 

(5.19) L = + 2 sech^ t. 

The latter’s analysis is greatly faciliated by the observation that 

{dt — tanht)((9t + tanhf) + 1 = L, while 
^ ^ {dt + tanht)(5t — tanhf) + 1 = d"), 

which reveals that for each integer n > 2 the kernel of L — is spanned by the functions 

(5.21) u±n{t) = {dt — tanht)e^"'* = (±n — tanht)e^"'*. 

The kernel of L itself is spanned by 

"Uoit) = — {dt — tanht)l = tanht and 

UQi{t) = {dt — tanht)t = 1 — ttanht; 
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(5.22) 



the Jacobi field uq is the normal projection of the Killing field generating vertical translations, 
while the Jacobi field uqi is the normal projection of the vector field generating dilations about the 
origin. Last we record the kernel of L — 1, though it plays no role in this construction, spanned by 

(5.23) u\{t) = {dt — tanht) sinht = secht and ui'(t) = {dt — tanht) = sinht + tsecht. 


which, multiplied by linear combinations cos 9 and sin0, respectively generate horizontal transla¬ 
tions and rotations about horizontal axes through the origin. 

From the above eigenfunctions of L, we see that for n > 2 the operator L — has a bounded 
Green’s function 


(5.24) 


Gn{t, s) 


1 j (n -|- tanh s) (n — tanh t) for t < s 

2n(l — n^) I — tanhs)(n + tanht) for t > s. 


The exceptional case n = 1 will not concern us here, but, for n 
Lu = f with trivial initial data, given by 


(5.25) 

with 



0, we will take the solution to 


(5.26) Gq {t, s) = tanh t — tanh s + (t — s) tanh t tanh s. 

Now suppose / G Cg^’"’'(]K). Then, setting 

= ^ fj 9) dO, 

^ 2 

(5.27) = — / f{t,9)cosn9d9, and 

Jo 

^ 2 

fn-{t) = - f{t,9)smn9d9, 

^ Jo 

the symmetries imply that fn-{t) for arbitrary n and /n,+ (t) for arbitrary odd n all vanish identi¬ 
cally. 

We then define 


(5.28) 


^ ^[joG^{t,s)%{s)ds ioi n = Q 

1 f-oo G 2 n{t, s)f 2 n,+ {s) ds otherwise 


It follows from 15.261 and 15.241 that 


(5.29) 


G 


supe \u 2 n{t)\ < 


teK 




0,0,7 


so in particular the series Yl’n=o ^2n(7) cos 2n9 converges in C^; we call the limit u(t). On the other 
hand, the series /o(t) + f‘in,+ it) cos2n9 converges distributionally to /. Accordingly u solves 

Gku = / weakly, but / G C;'|(“(]K), so in fact u G C'j^^“(K) is a classical solution to Cku = /. Thus, 
setting TZkI = u, interior Schauder estimates for Ck in conjunction with 
have constructed a right inverse 

(5.30) -Rk : C^’“’^(K) ^ Cc’“’^(K) 


5.29 


confirm that we 


to Ck- _ 

In turn, using also interior Schauder estimates for £;<;[*], we see that 

(5.31) 
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satisfies items (i) and (ii) of the proposition. 

The independence from ^ of and is already clear from 


5.6 


In proving the continuity in 


(^, we will emphasize the dependence of certain entities on the parameters by their inclusion, within 
brackets, in the notation. Of course Ck and TZk themselves are independent of the parameters, 
but motivated by 5.6 and recalling 5.4 we define 

(5.32) 
by 

(5.33) 


and we introduce 

(5.34) 

so that in particular 

(5.35) 


Ai[C] 

Aj[C](t, 0) = (sgn(t) arcoshr“^Tj(ri cosht), 0) 
^k[Q = A.[C]*“'>C,,A,[C]* : Cf’“g(K) ^ <’“g(K), 



(5.36) 


OO 

'^k[C\ = X] \^K {^K - ^k[C\^ 


n=0 


T^k-, 


whose convergence and boundedness, uniform in and whose continuity in C can be guaranteed 
by taking m large in terms of c, on account of the previous paragraph. Evidently, for each the 
map 

» — 1 . /^0,o,7 . ^2,a,-y / 


(5.37) 


K[crnK[c]K[cr 


. ^0,0,7 


(K) Ca’ ’^{K) 


is a bounded (by a constant which need not be assumed independent of Q right inverse for Ck, 
imposing trivial data at t = 0 for the rotationally invariant mode of the solution. On the other 
hand, the analysis of L above reveals that there can be just one such right inverse, namely TZk, 
because all higher modes annihilated by Ck grow too quickly in one direction or the other. Thus 
we hnd 


(5.38) 

concluding the proof. 


'PllT^icii][C,C]'Pc,^ = Ki[0,0]*-^TZK[C]i^i[0,0] 


-1. 


□ 


Toral solutions. Items (i) and (ii) of 4.33 imply 


(5.39) 

We set 

(5.40) 

and dehne C 

(5.41) 


lim nC^l 

vn-^OQ 


-1 


= A 


XT' 


C^’“’"(r[f]) = u|rw 


u G 


m' 


Qm . x ,— 


m 
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/^0,a,7 


im) 














by 


(5.42) 


£ 


r[i] 


= t: 




A'' T* 


To prepare for the next proposition we recall the definition 4.32| of wt, remembering in particular 
that wi and Wn are both identically 0, and we introduce Wi G for 1 < i < A^, defined by 


(5.43) 


Wi = Ip 


1 1 

lOim ’ 5£m 


op. 


Together these functions span the extended substitute kernel, modulo which we can invert TtIi] 
follows. 


Proposition 5.44. Given 0,7 G (0,1), c > 0, and positive integers k, i, and N, there exists 
mo > 0 —depending on the foregoing data—and there exists C > 0 —independent of c—such that 
for C,f,G [—c, I < i < N — 1, and m > mo, the map 

(5.45) ('r[i]c=c=o) ^ (r[l]c=€=o) 

has norm hounded by C, is constant in and depends continuously (with respect to the operator 
norm) on and moreover there exists a linear map 


(5.46) : { / G Cj“(r[i], x) | supp f C {p < 20£m}} ^ {T[i]) x M x M 

such that if f G C^’°^{T[i],x) is supported in {p < 20t'm} and {u,p,Jl) = TZY[i]f > then 

(i) Cr[i]U = f + pwi + JiWi; 

(ii) 11^112,a ,2 + \lA + \U\<C ||/||o,«; 

(Hi) the map 'iT2'R-T[i\Pc,i! projecting onto the component p of substitute kernel, is constant in (C,0/ 
(iv) the map 

(5.47) : { / G (r[i]c=«=o, x) | supp / C {p < 20£m}} ^ (r[l]c=^=o) 


has operator norm bounded uniformly in the parameters and m by C, is constant in and 
depends continuously (with respect to the operator norm) on (j; and 
(v) the map has norm bounded uniformly in the parameters and m by C, is constant 

in and depends continuously on (. 


Proof. Suppose / G Cg°‘{'T[i],x) has support in {p < 20im}. Extend its pullback T*f by 0, 
without renaming, to a function on the whole rectangle Tx,y, 0 ) and likewise extend Tfwi. Note 
that, via Ti, the symmetries of G require that Tf f extend to a function on invariant under 
reflections through the lattices of the construction and so doubly periodic, with periods 2X and 
2y. Observe that Tfwi enjoys these same symmetries, and in particular both functions satisfy 
periodic boundary conditions on Tx,y,o. Of course, the flat Laplacian on the torus obtained 
from Tx,y,o by identifying opposite edges has one-dimensional kernel consisting of the constant 
functions, so if we set 


(5.48) 


hx.Y.o'^i 


then T* (/ + pwi)\ is orthogonal to this kernel. The support of /, definitions of p and Wi, and 
dimensions of the rectangle guarantee \p\ < C ||/||(^o and 

(5-49) \\t: [p^ if + pw,)] < Cm ||/||^o , 
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so a solution uq G {'^x,Yfl-,9E) to 


(5.50) 
exists with 


(5.51) 


= T* [p'^ (/ + pwi)] 


llduol 


9e 


L^{9e) 


< Cm 


-1 


(JO , which in turn yields 


l^^ollco — 


-1 


llduol 


9E 


L^{9e) 


< c 


llduol 


9e 


< C 


CO 


L^(9e) 

Since / G C^’°‘{T[i],x)j we have uq G C'^’°‘{Tx yo, Oe), and in fact uq is harmonic outside 
{TtP < 20£m}. This region consists of one disc if T[i] is extreme and of two discs if T[i] is 
intermediate. Of course T^^T[i\ does not contain the center(s) of the disc(s), but the rectangle 
1’x,Y,o does, so we may define q to be the average value of uq on the centers (so simply the value 
at the single center in case i G {l,iV}). Then ui = uq — q is another solution of 5.50, since q is 
constant, and also satsifies 




co¬ 


if T[i] is extreme, then ui vanishes at the origin, but ui respects all the reflectional symmetries 
described above, so in particular is invariant under reflections through the origin, so being harmonic 
on the disc about it of radius 2 o 5 ii> accordingly decays there like the square of the distance to 
the center. In the intermediate case, however, we must perform an additional step to arrange the 
decay. Now, having subtracted q, the adjusted solution ui takes opposite values at the centers of 
the two discs in question. Define ~jl to be the value taken at the center where the function T*Vi 


(4.31) is —1. Then u = ui + pT*Vi vanishes at both centers and is still harmonic on both of the 
2 ^^-, ensuring the decay, but at the cost of altering the right-hand side of 


discs of radius 
that, recalling 4.32, u instead solves 


5.50 


so 


(5.52) 

and satisfies 

(5.53) 


^9E^ = Tt [p'^ if + pwi + pwi)] 


u 




-n-l 


4 , QE^ 2 

P 


<C\\f:C\T[i] 


Now the symmetries of T*f, T*Wi, T*Wi, T*p, and imply u = T- 
definitions xt = P^9e and = T*~^A 2 j.T* imply 


- 1 * 


u is ^-invariant, and the 


(5.54) 

so, using also 


5.53 


Cr[{\'^ = f + pwi + pwi, 

coupled with interior Schauder estimates for , we see that by defining TZ-pii 


according to TlT[i]f = u Y/e have established items (i) and (ii). 

The continuity and boundedness assertions concerning made at the beginning 

of the proposition, can be verified much like the analogous statements on catenoidal regions, by 
considering coordinate expressions for the conjugated operators and taking m large in terms of c. 
Item (hi) is trivial in view of the above construction of TZ-pii] the observation that 

( 5 . 55 ) T,[c,crr^^^w,or-^ 

acts as the identity on functions supported in {p < 20im}. Off this set each solution produced by 
is smooth and in fact harmonic, so in particular we get the estimate 

t 2' 


(5.56) 




m 


T[i\ n {p < 20^m},x, ^ 


< c 


0,a ’ 


but on the other hand, by taking m large in terms of c, we can ensure that 

,2' 


(5.57) 


V, 


-1 

C,? 
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is uniformly bounded and depends continuously on completing the proof. 


□ 


Global solutions. 


Proposition 5.58. Given a, 7 G (0,1) and initial surface data k, N, and parameter range c > 0, 
there exists C > 0 —independent of c—and there exists mo > 0 such that whenever C,? G [—c, 
and m > mo, there exists a linear map 




(5.59) 

such that if f £ = C = 0]) o,nd {u, {pi 

(i) C^u = + x; IJ'iWi + E 

(ii) V^lu + E 1^*1 + E ImJ < C I|/Ilo,o,.y; and 

(Hi) (v, 


t,N 


t,N-2 


m), (^ 2 , • • • , hN-i)) = HT'an 


^^u,{pi, - ■ ■ ,pn),{P 2 £" depends continuously on {(, . 


Proof. We will need the cutoff functions ^k, G C°°(S) defined by 

= ''P [20£m, 30£m] o p and 

(5.60) ^^^|lo.iE\(eU"7')'C|'l 

{if [rj“^ sech bi,T~^ sech 25j] o p on IC[i], 

and below we will casually identify these functions as well as the constant function 1 with the 
respective operators on spaces of functions that each defines by multiplcation, so that 1 acts as the 
identity, smoothly cuts off functions to have support contained in the union of the catenoidal 
regions, and 'b'r smoothly cuts off functions to have support contained in the union of the toral 
regions. We can also regard multiplication by these functions as extending to the entire initial 
surface functions originally defined on a single region. Additionally we will identify each region, 
catenoidal or toral, with its image under Q and will accordingly interpret TZ/cii], T^T[i]> 

Ct[i\ acting on these sets by the obvious extension enforcing the symmetries. 

Given / G set 


(5.61) 


then 

(5.62) 


N-l 

no = ('l'Kf)l/c[i] ! 

Z = 1 

A ^-1 

fl='^ ('^Kf)U[^] + (1 - f, 

i=l 

N 

ui = (/ilrw) > 

Pi = v^ 27 ^ 7 -[^] (/ilr[i]) for 1 < z < iV, and 
Pi = TTs'^rii] (/ilr[i]) for 2 < z < iV - 1; 

'^/= {uo + ui,{pi,--- ,pn),{p2,--- ^Pn-i)) > 


defines a map TZ : Cq’°‘’'^{Y,) 


^2, q ,7 


(E)x 


pN 


pN-2 


depending continuously on the parameters 


(via Pf^^) and bounded uniformly in them, according to the last two propositions, 5.15 and 5.44 
and to 15.91 
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Defining also L ; x x Cj"’^(S) by 


(5.63) 
we find 

(5.64) 


N N-1 

L{{u, (/XI,-- - ,/XAr),(/X2,--- ,/XAr_i)) = C^U - ^liWi - ^ /XjUJi, 


i=l 


i=2 


N N 

LTZf - f 4't] vri'^ 7 -[i] (/i|r[i]) + “ '^rw) (-/’ilrw) 


i=l 
N-1 

+ 

2=1 


2 = 1 


^ (/ix - '^^w) i^KfljCli]) - 
2=1 

The commutator factor in the hrst term is supported near the catenoidal boundary of the toral 


region, where by 5.44 the toral solution to /i has decayed rapidly, and of course the norm of fi is 


controlled by the norm of /. For the second and third term we can use the first two items of 4.33 
and of 4.16 respectively to bound the deviation of from its regional limits. We get 


(5.65) 


(^LTZ “ ^ C* (jn^ + m ^ + sech^ 6^ 


0 , 0:,7 ' 


where C is independent of m, c, the parameters themselves, and h. 

Thus by assuming b large in terms of C and then m large in terms of b and C, we can ensure 
LTZ is invertible on C'g“’"’^(S), with inverse continuous in the parameters (via T*^,^) and bounded 

uniformly in them by 2. Taking TZ = TZ concludes the proof. 

□ 

As an immediate corollary we obtain an estimate for the first correction to the surface as well as 
a statement of its continuous dependence on the parameters. 


Corollary 5.66. If 


(5.67) 

then the map 

(5.68) 


(tto, (/xi, - - - ,/XAr), (/X2, - - - - Y^ViW^ , 


(C,0 ^ (/^i, • • • ,hN),{h'2r-- ,hN-l)) 


is continuous, using the norm on the first factor of the target, and there exists C > 0 

depending on just k, i, and N such that for m sufficiently large in terms of k, i, N, and c 


(5.69) 




2,0,7 


+ ^ IT'-*! + ^ \hi\ < Cti- 


Proof. We note that T>i is manifestly continuous in the parameters and that the functions V^^H, 
and in fact vary smoothly on the initial surface and in the parameters. The continuity 

then follows from the proposition, as does the estimate, using also 4.43 (with a larger Holder 
exponent). □ 
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6. The main theorem 


We will need the following estimate for the nonlinear contribution 

( 6 . 1 ) Q[u]=n[u]-n[Q]-Cu 

that the perturbing function u makes to the mean curvature. 

Lemma 6.2. Given k, i, N, a,'y G (0,1), c > 0, and C > 0 independent of m and c, there 
exists mo > 0 sueh that whenever C, ? £ [^c, m > mo, and u : = ^ = 0] —)■ M satisfies 

||R|l 2 a 7 — C'ri, we have 


(6.3) 

Furthermore the map 




0,a,7 




(6.4) 


V-^Q[P-] : = e = 0]) X [-c,c] 


2N-2 


C°’^’^(E[C = e = 0]) 


V 


-1 




is continuous. 


Proof. If denotes the defining embedding of the initial surface with parameter values ((", ^), then 
V^^FL[P(^^^u] is simply the mean curvature of the perturbation by u of the embedding o P^^. 
Its value at p G = ^ = 0] thus depends smoothly on p, (C, 0; the 2-jet of u at p, implying 
the asserted continuity. 

Now hx p G = 0]. If we blow up the ambient metric gs to {p)9s and the 

perturbing function u to 


( 4 ’)'*^) tind 


(6.5) 




Jp 






where the latter two instances of Q indicate the nonlinear part of the mean curvature, relative to 
the superscripted ambient metric, of the perturbation by the bracketed function (in the normal 
direction corresponding to o) of the subscripted immersion. 

We can apply 4.24 to the blown-up situation to estimate the rightmost nonlinear terms appearing 


m 


6.5 on the ball Bp C S[(^ = ? = 0] with y-radius 1 and center p. Here the metric on Bp induced 


by the reparamerized immersion and rescaled ambient metric will be 


( 6 . 6 ) 


'PaPh 




where the subscripts on p and y emphasize that they are dehned on the initial surface with param¬ 
eter values Q and but then for j G {0,2}, the C^’°‘ norms on Bp induced by this metric and by 
yo,o will be equivalent, with ratio bounded above and below by positive constants independent of 
p, m, and c. 

Thus, noting that all derivatives of the second fundamental form of relative to the blown- 

up metric (and of course all derivatives of the curvature of the blown-up metric) are bounded by 
a constant independent of m and c, we apply the last two items of 4.24| to obtain from 6.5 the 
estimate 

(p)n:C2’“(Hp,yo,, 


(6.7) 




|QM : C'°-“(Hp,yo,o)|| < C (p) 

< Ce'^’^plQ{p)mfi"<P q‘P{p) ||r 


|2 

12 , 0,7 ' 
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where C may be different from the given one but is still independent of m and c. Therefore 

i+? 


( 6 . 8 ) 


p ^QMI 


0 , 0:,7 

for m large enough, as asserted. 




□ 


5.66 


We will need also a refinement of 3.21, after some preliminary definitions. Take (no,Ai,/i) as in 
and given v G C'g“’’^(S[C = i = 0]), let V = —TT'iTlp~‘^Q [uq + Vc^^^v]. Now define 


(6.9) Vi[v,C,,i]='^i + Pi + Vi 
for 2 < z < N — 1 , understanding Vi = Vn = 0, and define 

( 6 . 10 ) 


as the force in the direction K on the perturbation by uq + of fli (from 2.17). 


Lemma 6.11. Given a, 7 G (0, 1) and initial surface data k, I, and N, there exists c > 0 and there 
exists m-o > 0 such that whenever ^ > "lO; CLnd v G C|(S[C = ^ = 0]) satisfies 

II II ^ 


( 6 . 12 ) 


, we have, abbreviating Fi\v,Q,£f\ as 

Ti and Pfiv, C,^] as Pi, 



/ \ 


( P 1 + P 2 \ 


C 2 

kfw? 1 

Zj 


-2Lz-' 

P 2 + P 3 



2 'kt\ 


n 



Vc^-v 


\-7v-i — fF]^) 


\^Ar-l -1- Pn) 



< c, 


recalling Z : ^ —>■ ^ from 3.23, and 


(6.13) 


recalling H : 




/ Fi + Fn\ 


( \ 



k£m‘^ 

0 

2^-1 

P 2 








W/ 


V 0 j 


\Pn-i) 



< c, 


oTV-l 


oAf-l 


from 3 . 24 , Furthermore, for 1 < i < N, the maps (u,C7) 4Fi and 


are continuous, using the C'^’“(S,x) norm on the first factor of the domain. 

P roof. The conti nuity of Pi is clear from the definition of Pi and from the continuity statements 


m 


5.66 


6.2 


and 


5.58 


The force T) can be expressed as an integral along the image curves of 

-i„ 


dQi C S[C = ^ = 0] under the perturbation by + u of o P~^. Fixing a point p on such a 

curve, the image point and conormal there can be seen to depend smoothly on the p, the parameters, 
and the 1-jet at p of the perturbing function. Since the Killing field K and the ambient metric 
gs are themselves smooth, this establishes, in conjunction with the same references just cited, the 
asserted continuity of T). 


To check the estimates, first note that 


Pi — Pi 


< Cti, with C independent of m and the 


parameters. The proof will be completed by showing that 


iPi — fFi 


<Ct\ and then appealing to 


3.21 For the force comparison, set u = uq + Vc_,(^v, write Vu for the unit normal obtained from 
V by perturbing the initial surface by u, write \/\gfi for the area density of this perturbed 

surface, and write \/\g\ for the area density of the initial surface, as well as for the corresponding 
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length density on curves therein. Since K is Killing, we can use the first-variation formula and 
Green’s second identity to obtain 


— Ti = 


(6.14) 


/o. 


+ f - {u,K))n[u]y/\i\ 


+ / {i',K)Q[u]y/\g\ + / {u, K)r]uy/\i\ - / uri{u, K)y/\^\, 


in. 


IdUi 


IdUi 


recalling that rj is the outward conormal on 50^. 
From 


4.24 


(6.15) 


(i)) Il^lli 0 7 — estimates of ||A|| from 4.16 and 

V\9v] - VWl < Cm'^T^y/\g\, 


4.33 


we find 


so, using also the obvious boundedness of ||iF||^^, the estimate p ^ H^MHo < Cti, and the upper 
bound C{m? + |C|) for the x = P^9 area of flj, we can certainly bound the first term of 6.14 by 
Cm~‘^Ti, with C independent of m and the parameters, provided the former is chosen large in terms 
of the latter. Similarly, the difference between i^u and ly, parallely transported along the geodesics 
it generates, is controlled by ^1 -|- ||^||g IM^IIg and so in turn bounded by Cna^r'l, and of course 
|I1K|L_ is uniformly bounded, so we can bound the second term of 6.14 by Cm~‘^Ti as well. To 


19 s 


achieve the same bound for the third term we apply [6l2 


On the waist circles of the catenoidal regions we have \r]u\ < while the circles themselves 

have length of order ri, ensuring the estimate for the circular component (s) of the boundary in the 
On the rectangular component R, u satisfies periodic boundary conditions. 


fourth term of 6.14 


so j^rjUy/g = 0, but \{v,K) — 1| < Cm~‘^, |t/m| < Cmri, and R has length of order m~^, so we 
have Cm~‘^Ti as an upper bound for the entire fourth term. Finally, for the fifth term, on R we 
see \'q{v,K)\ < Cm~^Ti, while on the waist circle(s) we use the estimates 


(6.16) \rj{iy,K)\<\DK\^^ 

and |ri| < CrrRrl^'^, concluding the proof. 


+ 


\K\ 


9S — 


< Ct- 


-1 

1 


□ 


We can now prove the main theorem. 

Theorem 6.17. Fix 0,7 G (0,1). Given positive integers k, £, and N >2, there exist c,C,mQ > 
0 sueh that for every m > tuq there exist parameters 0 ,'^d a function u G 

such that ||rt|| 2 Q..), < Cti and the normal perturbation : S —)■ fey tt 

of the defining embedding t : S —?• is a minimal embedding. 

Proof. Set 


l 2,«,7 ^ T 


1"^"} X [-C,c]^ ^ X [-c,c]^ \ 


(6.18) il = |^;eCg2(S[C = ^ = 0 ],x) : 11^; 
and define J : B ^ B hy 


(6.19) 


c- 




k£mf 


^-1 


27rri 

k£w? 


Ti — Ti 


i+l 


7V-1 dvr _i /— ~ \ 

- ^ (T’i -|- Bi+i ) 

i=l Ti V 


N-1 
/ i=l 


Svr^r- 




where V and F abbreviate V[P(^^^v,C,,^] and C, ?] respectively. 
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By 5.66, 5.58 6.2, and 6.11 (which provides c), we are assured that J{B) C B as advertised. 
By the same references J is continuous, using the {^,x) norm on the first factor (and the 
usual ones on the second two). Moreover, B is in fact compact under this topology, because of 
the bound on elements of the first factor, and it is obviously convex. Thus Schauder’s fixed- 
point theorem ensures the existence of a fixed point vi to J. Set uo = —ttiTZ [p~^'H[0] — J2ViWi), 
vq = and u = uq + vo¬ 

lt follows that Ti = Vi = 0 for every i and that 


p — p ^?^[0] + Vy/ao + By^vo p + uq] + ^ ^ ptWi + 


'^PiWi 


( 6 . 20 ) 


= ^ ViWi + 
= ^ ViWi + 


E 


PiWi 


+ E 


PiWi 


PiWi 


but then Vi = pi Wi{i>u-,K), so, because the integrand has a sign, it is easy to see that the 

JUi 

vanishing force condition implies pi = 0 for every i. This sho ws minima lity, when ce follows the 
regularity since u G The estimate on u is clear from 5.66, 6.2, and 5.58, and in turn it 

implies the embeddedness. 

□ 
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